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We derive a microscopic equation of motion for the dynamical orientational correlators of molec- 
ular crystals. Our approach is based upon mode coupling theory. Compared to liquids we find four 
main differences: (i) the memory kernel contains Umklapp processes if the total momentum of two 
orientational modes is outside the first Brillouin zone, (ii) besides the static two-molecule orienta- 
tional correlators one also needs the static one- molecule orientational density as an input, where 
the latter is nontrivial due to the crystal's anisotropy, (iii) the static orientational current density 
correlator does contribute an anisotropic, inertia-independent part to the memory kernel, (iv) if the 
molecules are assumed to be fixed on a rigid lattice, the tensorial orientational correlators and the 
memory kernel have vanishing 1,1' = components, due to the absence of translational motion. The 
resulting mode coupling equations are solved for hard ellipsoids of revolution on a rigid sc-lattice. 
Using the static orientational correlators from Percus-Yevick theory we find an ideal glass transition 
generated due to precursors of orientational order which depend on Xo and <p, the aspect ratio and 
packing fraction of the ellipsoids. The glass formation of oblate ellipsoids is enhanced compared to 
that for prolate ones. For oblate ellipsoids with Xo < 0.7 and prolate ellipsoids with Xo > 4, the 
critical diagonal nonergodicity parameters in reciprocal space exhibit more or less sharp maxima 
at the zone center with very small values elsewhere, while for prolate ellipsoids with 2 < Xo ;$ 2.5 
we have maxima at the zone edge. The off-diagonal nonergodicity parameters are not restricted to 
positive values and show similar behavior. For 0.7 < Xo < 2, no glass transition is found because 
of too small static orientational correlators. In the glass phase, the nonergodicity parameters show 
a much more pronounced q-dependence. 

PACS numbers: 61.43.-j, 64.70.Pf, 63.90.+t 



I. INTRODUCTION 

The experimental and theoretical investigation of sys- 
tems with self-generated disorder has mainly been de- 
voted to simple and molecular liquids P, |J. In their 
supercooled state at low temperatures or high densities, 
liquids exhibit nontrivial dynamics, often called glassy 
dynamics. Decreasing the temperature T or increasing 
the number density n may result in a glass transition. 
The physical origin of glassy dynamics and the glass tran- 
sition is the formation of a cage by the particles. For 
not too low temperatures and not too high densities, the 
cage's lifetime is finite, i.e. a particle can escape with a fi- 
nite probability. However, if the lifetime diverges, e.g. at 
a critical temperature, the particles remain localized in 
their cages. In that case an ideal glass transition occurs 
at that temperature. 

There is no general agreement about the theoretical de- 
scription of the glass transition. From a practical point of 
view one often considers the so-called calorimetric glass 
transition temperature T g as the temperature at which 
a supercooled liquid becomes a structural glass. But T g 
depends on the cooling process. Therefore, it is not well 
defined. Besides T g , there exist two better defined char- 
acteristic temperatures, which are T c , the mode coupling 
transition temperature, and Tjc, the Kauzmann temper- 
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ature. At Tk the excess entropy of a supercooled liquid 
with respect to its crystalline phase disappears. This is 
a very old concept which only recently has been put onto 
a microscopic basis by the replica theory for structural 
glasses (see |2| and references therein). The existence of 
a purely dynamical glass transition at a critical temper- 
ature T c has been suggested about two decades ago 
This approach is based on mode coupling theory (MCT). 
MCT describes the cage effect (as explained above) in a 
self-consistent way. At T c a transition from an ergodic to 
a nonergodic phase occurs. Close to T c the relaxational 
dynamics of e.g. the density fluctuations, exhibits two 
time scaling laws with relaxation times which diverge at 
T c . For more details and comparison of the MCT pre- 
dictions with experimental and simulational results the 
reader is referred to Refs. 0, IE El El ■ A review of 
MCT, replica theory of structural glasses and a selection 
of phenomenological theories is given in Ref. 0] • 

Glassy behavior of systems with self-generated disorder 
is not restricted to liquids. There exist so-called molecu- 
lar crystals 0] were molecules are located at sites of a pe- 
riodic lattice. At higher temperatures their orientational 
degrees of freedom (ODOF) may be dynamically disor- 
dered, i.e. ergodic. This phase is called the plastically 
crystalline phase 0] . Decreasing temperature lowers the 
lattice constants which in turn leads to an increase of the 
steric hindrance between the ODOF. This may result in 
the formation of an "orientational" cage in which the ori- 
entation of a molecule is captured on a certain time scale, 
quite similar to liquids. If this time scale diverges the 
plastic crystal undergoes an ideal orientational glass tran- 
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sition. The corresponding phase is called glassy crystal. 
That such a glass transition really occurs has been proven 
experimentally several decades ago. First systematic ex- 
perimental indication for the formation of glassy crystals 
has been given in 1974 for several molecular crystals [l3[ . 
Since then a lot of glassy crystals were found. With- 
out claiming completeness the most intensively studied 
molecular crystals are cvanoadamanta ne 1 1 4 Il5l ll6L Il7| , 
chloradamantane ^3 and ethanol [T^ . |20|. Ethanol has 
the big advantage that it can form either a supercooled 
liquid, a structural glass, a plastic crystal, a glassy crystal 
and an orientationally ordered crystal within a small tem- 
perature interval around 100 K. Therefore, it has been 
investigated experimentally to explore the role of transla- 
tional degrees of freedom (TDOF) and ODOF for glassy 
behavior |20j. These experiments have shown that the 
ODOF of molecular crystals exhibit quite similar glassy 
behavior than conventional supercooled liquids. Addi- 
tionally, comparing different molecular crystals with each 
other, similar glassy behavior was found 15| . These sim- 
ilarities also include dynamical heterogeneities (2l| . The 
largest deviations of molecular crystals from supercooled 
liquids were observed in dielectric spectroscopy. The for- 
mer exhibit a rather weak excess wing, or even no such 
wing, in contrast to supercooled liquids [T^. 

An interesting model for molecular crystals has been 
studied some years ago. The molecules were approxi- 
mated by infinitely thin hard rods with length L which 
were either fixed with their centers on a fcc-lattice j2^ 
or with their endpoints on a sc-lattice |2^| . MD- and 
MC-simulations (23, respectively, have shown the exis- 
tence of glasslike dynamics. Particularly a critical length 
l c — L/a (a is the lattice constant) has been determined 
at which an orientational glass transition occurs [22I |23| . 
However, this transition is not sharp, in close analogy to 
supercooled liquids. The system of infinitely thin hard 
rods is particularly interesting since there are no static 
orientational correlations. Consequently, glassy behavior 
does not originate from growing static correlations, but 
results from entanglement which leads to a "dynamical 
cage" . 

As far as we know there is no microscopic theory which 
describes glassy behavior of molecular crystals with self- 
generated disorder. For mixed crystals [24j. i.e. crys- 
tals with quenched disorder, a microscopic theory for the 
glass transition has been worked out |25|. This theory 
is based on MCT and takes into account ODOF and 
TDOF, i.e. lattice displacements, as well as translation- 
rotation pt| coupling. The displacements are crucial, 
since the statistical substitution, of e.g. CN molecules 
in KCN by Br atoms, leading to the well-known mixed 
crystal compound (KBr)i_ K (KCN)^ |24|, generates ran- 
dom displacements. Due to the translation-rotation cou- 
pling, these random displacements induce random fields 
acting on the ODOF. MCT was also applied to spin 
glass models, where the coupling constants between spins 
are at random [^3 ■ Unfortunately, the quality of MCT 
predictions for mixed crystals and spin glasses has not 



really been tested, in contrast to supercooled liquids 

Since MCT has been very successful HHQ,|lli to de- 
scribe glassy dynamics of supercooled liquids, and since it 
has also been applied to mixed crystals and spin glasses, 
it is natural to derive MCT equations for plastic crystals, 
as well. This will be done in Sec. [HI The calculation of 
the glass transition line and the critical nonergodicity 
parameters from the MCT equations will be presented 
in Sec. IIIII for hard ellipsoids on a sc-lattice. The final 
section contains a discussion of the results and some con- 
clusions. More technical details are put into four appen- 
dices. 



II. THEORETICAL FRAMEWORK 

In this section we will describe how MCT equations 
can be derived for molecular crystals. The strategy is 
quite similar to that for simple liq uids 0, |(|, molec- 
ular liquids of linear molecules |28l l29l | and arbitrary 
molecules |30J . The introduction of the microscopic ori- 
entational density, the corresponding current density and 
their time-dependent correlators will be described in sub- 
section lTl Al Then, in subsection lll Bl wc apply the Mori- 
Zwanzig projection formalism |3lll32| to derive an equa- 
tion of motion for the time-dependent orientational cor- 
relators. Following MCT for liquids, the memory kernel 
is then approximated by a bilinear superposition of the 
time-dependent orientational correlators. 



Microscopic orientational densities and their 
correlators 



We consider a Bravais lattice with N lattice sites. 
Since the experimental and simulational results for super- 
cooled molecular crystals O EE EE E3, E3, E3, 123l 
have demonstrated that glassy behavior can occur due 
to steric hindrance of the ODOF, we restrict ourselves 
to a rigid lattice, i.e. we neglect the translation-rotation 
coupling. Then, the increase of steric hindrance by de- 
creasing temperature can be accounted for by either a 
variation of the lattice constant or equivalently by an in- 
crease of the size of the molecules. At each lattice site 
we fix a molecule. The natural way is to fix its center 
of mass. All molecules are assumed to be identical and 
rigid, as well. We will consider linear molecules only. 
Generalization to arbitrary molecules can be performed 
like for molecular liquids |30| . 

It is also obvious that the ODOF are best described 
in a molecule-fixed frame with its origin coinciding with 
the lattice site. In principle one could choose any other 
reference point [2l| . But this would artificially introduce 
TDOF, besides the ODOF. Using the former choice, the 
ODOF of the n-th molecule at the site with lattice vec- 
tor x n is given by the angles tt n = (0 n ,(f) n ). The third 
angle Xn with respect to the symmetry axis is irrelevant 



3 



for the glassy dynamics. The moment of inertia for the 
axes perpendicular to the symmetry axis is denoted by 
/. The interaction between the molecules is given by 
V(f2i, . . . ,f2jv) and the classical dynamics follows from 
the classical Hamiltonian 



H({9 n ,<j>n},{Pe n iP^ n }) 



1 N 

- y 

21 ^ 



n=l 



Pi 



Pi 



+ v({e n ,M), (i) 



where pe n and are the momenta conjugate to n and 
cj> n , respectively. 

Next we introduce the microscopic, local orientational 
density 



p n (Sl,t)=6{Sl\n n (t)) : 



(2) 



with d(Q\Q') = (sm6)- 1 S(9 - 6')5((j) - </>') and Q n (t) 
the classical trajectory of the n-th molecule. The one- 
molecule orientational density p^(Q) is given by 



pM(ft) - (<5(Q|fi n (t))) , 



(3) 



which is independent on t and, for identical molecules, 
also on n. ((•)) denotes canonical averaging with respect 
to initial conditions in the 4iV-dimensional phase space. 

Taking the time derivative of p n (f2, t) leads to the con- 
tinuity equation. 



(4) 



Here, Ln^ is the angular momentum operator acting on 
fl n , and 

j n (n,t)=(* n (t)6(n\fi n (t)) = (* n (t) Pn (si,t) (5) 

is the corresponding orientational current density, which 
involves the angular velocity w n (t). We also introduce 
the "longitudinal" orientational current density 



j n (Q,t) = Lq„ • j„(Q,t) 



(6) 



With these quantities we can define the time-dependent 
orientational correlators: 

GW(n,fi',t) = (8p n (n,t)s Pnl (n')) (7) 

of the local orientational density fluctuations 

s Pn (Q,t) = Pn (n,t)-(p n (n,t)) = Pn (n,t) - , 

(8) 

as well as 

j nn , (n, n', t) = 4tt ( j„(n, t) j n . (o') ) . (9) 



Here, Sp n (Q) = Sp n (Q,0) and j n (Q) = j n (fl,0). 

Similar to molecular liquids |28l l29l |. we expand the 
orientation-dependent functions with respect to spher- 
ical harmonics Yx(fl), A = (Irn), as already done for 



the static correlators |33j[51|. This allows to repre- 
sent any functions f n (Sl, t) and F nn i (f2, 0', t) by their A- 
and Fourier-transforms. The corresponding transform of 
G nn i{Q,,Q! \t) leads to the intermediate scattering func- 
tions 

Sxx'(q,t) = — (6p* x (q,t)5p x ,(q)) (10) 
and the corresponding current density correlators 

Jxx-{^t) = —{3lMh' '(q)>. (ii) 

These correlators form matrices S(q, t) — (Sw(q,t)), 
etc. The wave vectors q are restricted to the 1. Brillouin 
zone (l.BZ), due to the lattice translational invariance. 
The correlators Sw(q,t) form a complete set. For ex- 
ample the neutron scattering function S , ne utron(q, t) can 
be expressed by {5AA'(c[,t)} using the scattering lengths 
of the molecular sites [36j |. 

Note that the correlators (|1U|I and Ijllfl vanish for I = 
and/or /' = 0, due to the absence of TDOF. The symme- 
tries of the orientational correlators discussed in Ref. j33| 
also hold for the time dependent quantities. They will be 
applied to reduce the number of independent correlators. 



B. Mode coupling theory 

The goal of this subsection is to derive an approxi- 
mate equation of motion for the intermediate scattering 
functions Sw (q, f) of molecular crystals. Due to the 
rigid lattice, only ODOF are involved. In case that the 
steric hindrance is large enough the orientational den- 
sity fluctuations Sp\(q,t) contain slow parts. Choosing 
6px(<l,t) and the "longitudinal" current density jx(<l,t) 
as slow variables, we can apply the Mori-Zwanzig formal- 
ism [3lTl3^ to derive an equation of motion for S(q, t): 



S(q,t)+JS- 1 (q)S(q,t)- 



dt' M(q, t — t') J S(q, t') = . (12) 



The notation _1 means the inverse of the 1,1' > block 
of the respective matrix, i.e. the inverse with respect to 
the subspace of non-constant functions in angular space. 
This is because the first rows and columns of these ma- 
trices vanish. The only exception from this rule is d" 1 
in App. |D] The prefactor J S _1 (q) in Eq. is related 
to 



^(q) = S- 1 / 2 (q)JS-^(q), 



1/2/ 



(13) 



which is the square of the symmetric microscopic fre- 
quency matrix (f2^^/(q)). It depends on the static orien- 
tational correlators S(q) and on J e (</\A'(q)), which is 
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independent of q (see App. [XJ. The matrix elements of 
the memory kernel M(q, t) are given by 

Air 

M xx ,(q,t) = — ((£jx(q))*\Qe- iQCQt Q\ £j x ,(q) ) , 

(14) 

the correlations of the fluctuating forces Q | £j\(q) ) . £ 
is the Liouville operator and Q = 1 — P p — Pj (see Eq. 
([Bill)) projects perpendicular to the slow variables 5p\ (q) 
and j x {q). 

In a final step we perform the mode coupling approx- 
imation for the slow part m(q, t) of J -1 M(q, t) J -1 , 
which enters Eq. J22J, yielding (see Appendices IBIDfl 



"»A*(q,t)«2^E E' E' ( 15 ) 

^(qAA'|qi AiAi; q2A 2 A' 2 ) S Al a; (qi, *) SasA!, (<J2, *) • 
The vertices are 
F(qAA'|q 1 A 1 A' 1 ;q 2 A 2 A' 2 ) = 



(4tt) 2 E ( J \ 



AA :1 



A3 A 



3/\ 3 



X!' w(qA 3 |qiAi;q 2 A 2 ; A") 



x 



E' <qA^|qiA' i; q 2 A 2 ;A'") 



L A" 
where 



(J" 



' AiA' 



(16) 



w(qA|qiAi;q 2 A 2 ; A") = 

= 6(i"i 2 C(l"hl, m"m 2 m) c x „ Xl (q x ) + (1 «- 2) , (17) 



(2i + l)(2/ / + l) 
21" + 1 



(18) 



1 + (-1) 1 \/7(7TT) y/i"(i" + 1) c(K'i", 101) , 



and 



■Aa' 



Jaa' 



(19) 



is the inertia and temperature-independent part of Jw . 
5^ ' denotes summation such that qi + q 2 = q + Q, 
with Q a reciprocal lattice vector, and indicates 
summation over all A 7^ (00). C(ll'l" ,mm'm") are the 
Clebsch-Gordon coefficients and caa' (q) the direct corre- 
lation function matrix elements. 

The result, Eqs. I|15l) - (|19l) . has a striking similarity to 
the slow rotational part rn^f, (q, t) of the memory ker- 
nel for molecular liquids 28, 29]. This is not surprising. 
Particularly u(qA|qiAi; q 2 A 2 ; A") are identical, up to a 



factor + 1). This similarity originates from the fac- 
torization of a static three-point correlator described in 
Appendix El It is this approximation which leads to the 
rather simple result, Eqs. (|16fl and 117|) . for the vertices. 
Of course, taking the static three-point correlator from a 
simulation would make this factorization approximation 
unnecessary. However, we do not expect any qualitative 
influence using our approximation instead of the correct 
simulational result. Such an influence is only to be ex- 
pected for systems which interact through three-, four- 
. . . body potentials. Since the molecules are fixed on lat- 
tice sites, such covalent bonds may be less important for 
molecular crystals. 

There are four main differences with respect to MCT 
for molecular liquids. First, the tensorial MCT equations 
for molecular liquids are first order integro-differcntial 
equations which can not be transformed to a second order 
integro-differential equation, like Eq. i|12fl . Second, for 
molecular crystals toaa' (q, t) contains a sum over recip- 
rocal lattice vectors Q. Therefore, the sum over qi , q 2 in- 
volves Umklapp processes. Third, due to the rigid lattice, 
only the 1,1' > matrix elements are nonzero. Fourth, 
the static current density correlator Jaa' in Eq. (|T5|l docs 
not cancel completely, as it does for a molecular liquid. 
There remains an anisotropic part Jaa', which equals 
1(1 + 1)S\\> for a liquid and is defined in Eq. JT^J. Jaa', 
can be related to the A-transform p { ^ } of the one-molecule 
orientational density pW(fi), which is needed as an input 
(for details see App. |A"|) . 

There is no explicit dependence of the kernel m(q, t) 
on T and /. On a large time scale we can neglect S(q, t) 
in Eq. IJ12H . The remaining equation does not involve any 
inertia effect, i.e. the glassy dynamics depends not on /, 
except for fixing the time scale. 

The vertices V^qAA'lqiAiA^; q 2 A 2 A 2 ) depend on J and 
the direct correlation function c(q) only, where c(q) is re- 
lated to the static orientational correlators S(q) by the 
Ornstein-Zernike (OZ) equation for molecular crystals 



S(q) = 



D 1 7-c(q 



4tt 



(20) 



This equation is similar to that for molecular liquids pa . 
|33[53, with exception of the appearance of D. D\\' 
is the A-transform of D(Q,Q') = 4tt S(fl\Q') - 

ff2) p^ (CI') ] and can be expressed by the A-transform 

of pW(n), too. 

This discussion makes clear that the closed set of MCT 
equations l|12|) - (|19|l requires two different static input 
quantities, the one-molecule quantities D\\r,J\\i which 
can be expressed by {p^} and the two-molecule correla- 
tors S\\> (q) , or ca A' (q) ■ Note also that we have neglected 
contributions to m\\i (q, t) coming from the fast part of 
Sp(q,t) which leads to a damping term in Eq. I|12|) . On 
a long time scale this has no influence. 

Here we restrict ourselves to the investigation of the 
orientational glass transition itself. For this we introduce 
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the nonergodicity parameters (NEP, not normalized!) 

F AV (q) = lim Sxx'fat). (21) 

t — *oo 

In the limit t — ► oo, Eq. Ijl2l) leads to 

S-^FCq) [S(q)-F(q)]- 1 =JP[F(q)] (22) 

where 

J"[F(q)] = lim m(q,i) . (23) 
Eqs. J121) and fflft 

are matrix equations for 1,1' > 0, 
since the first columns and rows of the involved matrices 
vanish, due to the rigid lattice. 

III. RESULTS FOR HARD ELLIPSOIDS 

After having derived equations of motion for the ori- 
cntational correlators S\y (q, t) , their time dependence 
could be calculated numerically. Although the mathe- 
matical structure of the MCT-equations (fl2*|l . dHJ-dS) 
is identical to that of more-component liquids, the nu- 
merical solution is hampered due to the anisotropy of the 
lattice, in contrast to liquids. Because of this anisotropy 
the correlators also depend on the direction of q. For 
liquids it has turned out that the restriction to several 
hundreds of values for q = |q| leads to rather precise so- 
lutions of the MCT-equations. We will see below that we 
have to choose several thousands of q vectors within the 
first Brillouin zone. In addition, in comparison to molec- 
ular liquids, there are more independent correlators for 
each pair (/, I'), increasing the number of equations even 
more. Consequently, a numerical solution will require ei- 
ther an improvement of the numerical code |38j usually 
used for the numerical solution of the MCT-equations 
and/or further simplification of S\y(q), e.g. neglect- 
ing the dependence on q/|q|. Since we want to avoid 
such type of approximations we will restrict ourselves to 
the calculation of the glass transition point and the cor- 
responding critical nonergodicity parameters and leave 
the solution of the time-dependent equations for future. 
Nevertheless, their identical structure to that of liquids 
already ensures the validity of, e.g., the two time scaling 
laws pj for molecular crystals as well. 

In order to solve Eq. I)22ll we have chosen hard ellip- 
soids of revolution, fixed with their centers on a sc-lattice 
with lattice constant equal to one. The symmetry axis 
of the ellipsoids has length a and the length of the per- 
pendicular axes is b. Replacing linear, rigid molecules by 
hard ellipsoids is probably not a bad approximation since 
the steric hindrance is qualitatively the same. In addi- 
tion, the choice of hard ellipsoids has two advantages: 
First, we have already calculated the static orientational 
correlators Sw(q) for 1,1' < 4 within the Percus-Yevick 

(PY) approximation, and (which yields Jaa' and 
D\\/ ) for / < 8 by MC-simulations. Second, we have 



recently solved the MCT-equations for a molecular liq- 
uid of hard ellipsoids [3(j which allows to compare the 
conditions for the appearance of the ideal glass transition 
for the ellipsoids on a lattice and in their liquid phase. 
This comparison will allow to estimate the qualitative or 
quantitative role of TDOF for the freezing of the ODOF. 
Furthermore, the ellipsoids' head-tail symmetry leads to 
a decomposition of Eq. i|12|) and therefore of Eq. 1221) into 
a closed set of equations for Sw (q, t) and F\\> (q) , re- 
spectively, for I, I' both even and a set of these quantities 
for 1,1' both odd. All correlators with / even and V odd 
or vice versa are zero. The set of equations for 1,1' both 
even is closed because the memory kernel only contains 
correlators with I and V even. This is in contrast to the 
equations for /, /' odd. The corresponding memory kernel 
contains a bilinear coupling of correlators with 1,1' even 
with those where 1,1' are odd. It is easy to prove that 

4 s A l(i)^S A v(q,i), (24) 

for 1,1' both odd. The "self" correlator S x %(t) is the A- 
transform of G„ n (f2, fi', t), up to a prefactor. In contrast 
to this the "self" correlator with I, I' even is given by 

S AX'(*) = ^ E SAA<(q,i). (25) 

qGl.BZ 

Similar relations hold for the NEP. 

There is a technical disadvantage connected with the 
hard body potential between the ellipsoids. The inver- 
sion of the static correlators, occurring e.g. in the projec- 
tors (cf. Ea.Q14[ll and vertices (cf. Eq. (|16fl ). needs some 
caution. Readers interested in this point are referred to 
Refs. l33[ [4(1 . We stress that the second paper cited 
in Ref. |33j|and |40| contain additional technical infor- 
mation, particularly the discussion of those mathemati- 
cal problems related to hard core interactions. However, 
these details will not be needed in the following. 

The numerical solution of Eqs. (|22|l . i)23|) requires a 
restriction of the matrices to /, I' < i max - We have 
chosen ^ max = 4. The q-vectors are discretized, i.e. 
for the a-component of q we have chosen q a = v a ^, 
v a = -f + I,-- - ,(),•■• ,f - 1 with M = 32, 
which makes a total of 32768 q-vectors. Due to the 
point symmetry of the lattice, the number of indepen- 
dent FAA'(q) can be reduced. For more details the 
reader is referred to Refs. [H, fipj . The solution of 
Eqs. i|22|) yields the NEP F\y(q) and the corresponding 
normalized quantities f\y(ci), where the normalization 
/AA'(q) = ^AA'(q) [S , AA(q)5 , A'A'(q)] _1/2 has been cho- 
sen. Varying the aspect ratio Xq — | and the volume 
fraction ip = ^ab 2 , we have located the glass transition 
line, at which for I, I' both even a nontrivial solution for 
F\\'(q) bifurcates. 

This has been done by approaching the glass transition 
point from the glass side. Fig.^shows an example where 
^2i,2i(q = 0) is represented as function of a for fixed 
b. The fit with a square root, predicted by MCT for a 
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1857 5.858' 5.859 5.86 



a 

FIG. 1: a-dependence of the NEP F^i^i (q = 0) in the vicinity 
of the glass transition point for fixed b = 0.24 and different 
values of a (squares). Also shown is the square root fit 47 + 
400\/a — 5.85688, leading to the critical values (a c , -£21,21(0)) 
(black dot). Note the rather large prefactor of the square 
root. 



type B-transition [5J allows to locate the glass transition 
point a c (b) for fixed b = 0.24 up to a relative deviation 
better than 1CT 4 ! The NEP F 2 i, 2 i(0) for a = 5.857 (see 
Fig. I1U|) deviates less than ten percent from the critical 
NEP F 2 C 121 (0) = 47 at a c = 5.85688, and no qualitative 
change is to be expected on further approach towards a c . 



A. Phase diagram 

The glass transition line ip c (Xo) obtained in this way 
is shown in Fig. [21 Fig. [21 also contains the equilibrium 
phase transition line ip cq (Xo) from MC-simulations, and 
the line ^py(^o)- Also shown is the curve (^extra(^o) 
where the extrapolated OZ/PY static orientational cor- 
relators for Xq > 4 at the zone center diverge. ip eq (Xo) 
and ippy(Xo) were obtained from the corresponding lines 
a cq (b) and apy(6) of Ref. [33|. Finally, the solid line with 
the cusp at Xq = 1 is the location of all <p(Xq) at which 
the rotators start to interact. 

At </? eq (Ao), an equilibrium phase transition from a 
(dynamically) disordered to an orientationally ordered 
phase occurs. The line ypY^o) locates the (Xq, y)-pairs 
for which the iterative numerical procedure to solve the 
OZ/PY equations becomes unstable. This is associated 
with some of the maxima of S(q) becoming very large, 
giving evidence of a divergency. Since the Y -dependence 
of ipeq and tppY is qualitatively similar, this behavior may 
indicate an equilibrium phase transition, as speculated 
for a liquid of hard ellipsoids |4l|. However, in contrast 
to the latter, the deviation of ifpy(Xo) from <£> eq (Yo) is 
much larger, especially for prolate ellipsoids. Using the 
static correlators from OZ/PY-theory as an input for the 
calculation of the NEP from Eqs. lt2"2"|) and (|2"5|) we have 
found a glass transition for Xq < 0.7 and 2 < Xq < 2.5, 
only. For 0.7 < Xq < 2 and Xq > 4 the system is ergodic 
for all ip < (Ppy(Xq), because the static correlators at 
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FIG. 2: Phase diagram of hard ellipsoids on a sc lattice. 
Shown is the curve below which the ellipsoids are free ro- 
tators (full line with cusp at Xq = 1), the equilibrium phase 
transition line ^eq(-^o) from MC simulations (full thin lines), 
the line <ppy(Xo) of highest densities to be reached by numer- 
ical solution of the OZ/PY equations (dash-dotted lines), the 
line </'extra(Yo) for Xo > 4, where the extrapolated OZ/PY 
results diverge (dotted line) and the MCT glass transition line 
ip c {X ) (full thick lines). For X > 8, (p e xtra(^o) and tp c (X ) 
are almost identical. • denote the five state points for which 
the nonergodicity parameters presented in Figs. 131121 were 
calculated. 



the Brillouin zone center and/or edge are too small. Un- 
fortunately, the iterative procedure to solve the OZ/PY 
equations for 0.7 < Xq < 2 and Xq > 4 becomes un- 
stable for ip > (Ppy(Xq). Therefore, we have decided to 
extrapolate the static correlators to ip > (Ppy(Xq). This 
extrapolation is guided by the physical assumption that 
long range orientational order should occur at the line 
<y3cxtra(Yo). It only works for Xq > 4, but not for the 
gap in between X 0.7 and X ss 2. Accordingly, the 
missing glass transition line for 0.7 < Xq < 2 first of all is 
based upon the lack of the static input. Since the ergodic 
and nonergodic phase are separated by a critical line of 
Type-B transitions 0,0 j Pc{Xq) can not terminate at 
Xq w 0.7 or Xo ~ 2. There exist two possible scenarios 
for (p c (Xq) within this gap. First, (p c (Xo) converges to 
tp c (Xo — 1) from above and below Xq — 1, with a pos- 
sible cusp at Xq — 1. Second, <p c {Xo) — > v 3 max(Xo t ) f° r 
Xq -> Xq with 0.7 < Xq < 1 and 1 < Y+ < 2, where 
Pm&xiXo) is the maximum possible volume fraction of 
an orientationally disordered configuration for given Xq. 
The second scenario would imply that there is no glass 
transition for Xq < Xq < Xq , i.e. for ellipsoids which 
are not sufficiently aspherical. 

The non-monotonous behavior of ippy{Xq) for pro- 
late ellipsoids with 1 < X < 4, which induces a non- 
monotonicity of </j c (A"o), seems to be an artefact of the 
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PY approximation, as our MC results for hard prolate 
ellipsoids suggest, though the static orientational corre- 
lators from OZ/PY theory are qualitatively correct, any- 
way [33|. 

If it is true that the divergence of the PY solutions cor- 
responds to an equilibrium phase transition, this implies 
that the ideal glass transition is driven by the growth 
of some S\\> (q) at the zone center or/and edge due to 
the growth of the orientational order, as will be seen in 
the following figures. This is quite similar to the cen- 
tral peak phenomenon above the equilibrium transition 
temperature at structural phase transitions of first and 
second order • The central peak can be interpreted as 
a quasi-nonergodic behavior and has also been described 
by MCT |H. 

The freezing of the 1,1' odd correlators occurs beyond 
the 1,1' even glass transition line and is treated in sub- 
section [^^Ql 



B. Critical nonergodicity parameters 

The critical NEP F xx ,(q) and the normalized critical 
NEP f XX r{q) together with the static orientational cor- 
relators are shown in Figs. EH El El EH and for oblate 
and prolate ellipsoids, respectively, along the three highly 
symmetric directions in reciprocal space from the zone 
center to its edge. For each of the three directions and 
each matrix element, a separate subfigure is provided, 
where the indices Iml'm' are displayed at the top of each 
figure. We have restricted our illustrations to the diago- 
nal elements I = I' = 2, m = m' = 0, 1, 2 and I = I' = 4, 
to = m' = 0,1,2,3,4, and the off-diagonal elements 
I = 2, I' = 4, m = m' = 0,1,2. By the symmetries 
of the cubic lattice, these correlators are all real. The 
scales on the l.h.s. of each tableau belong to S\y (q) and 
F xx , (q), those on the r.h.s. to f xx ,(q). Note the different 
scales of the axes for different values of to = to'. 

For two pairs (a, b) we also present the corresponding 
tensorial quantities in real space. Figs. 0] and 00 show 
log-lin representations of the direct space static orienta- 
tional correlators G xyz> w and the corresponding NEP 
F x V z,\\> = finij—Kx, G X yz,w(t) along lattice directions 
of high symmetry, i.e. xyz = OOn, Onn and nnn for 
n = 0, 1, . . . , 8. Along these directions, all G xyZt \\> and 
F£ xx , are real, too, for AA' as above. Note that a step 
An = 1 corresponds to different lengths in direct space, 
namely 1, y/2 and V3 for the different lattice directions. 
For each to = to' and each lattice direction, a separate 
figure is provided and a logarithmic plotting has been 
chosen for positive and negative values of G xyZt \y and 
F£ xx , separately, i.e. the negative values are presented 
as - In \G xyz .\\>\ and - In \F£ yzXX ,\, respectively. The 
values of xyz, Iml'm' are included in each subfigure. 

Figs. EE present the NEP for I = I' = 2 and I' = 2,4 
for oblate ellipsoids with a = 0.08 and b — 1.412, which 
yields (X ,<p) = (0.0567,0.0835). In comparison to liq- 
uids, the NEP possess less structure in g-space. For 
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FIG. 3: q-dependence of the nonergodicity parameters 
F XX /(ci) (dotted lines), the normalized ones /xv(l) = 
^AA'(q)[SAA(q)Svv(q)r 1/2 (thick grey lines) and of the 
static structure factors S XX i(q) (solid lines) for I = I' = 2 
and m = m' = 0, 1, 2. Results are shown within the first Bril- 
louin zone along the three highly symmetric reciprocal space 
directions for oblate ellipsoids with a — 0.08 and b — 1.412, 
i.e. (X ,(p) = (0.0567,0.0835). 
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(solid circles) and the static orientational density correla- 
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along the three highly symmetric direct lattice directions 
for oblate ellipsoids with a = 0.08 and b = 1.412, i.e. 
(Xo,ip) — (0.0567,0.0835). Shown are the diagonal corre- 
lators for I = I' = 2, m = m' = 0, 1, 2 and xyz = OOn, Onn or 
nnn for n = 0, 1, . . . , 8. 
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FIG. 5: n-dependence of the ratio -Foon,2i.2i/C?oon,2i,2i of 
the NEP Fobn 21,21 an d his static counterpart Goon,2i,2i for 
oblate ellipsoids with a = 0.08 and b = 1.412, i.e. (Xo,<p) — 
(0.0567,0.0835). 



responding G xyZt \\', as can be seen from G nnn , 20,20 and 
Fnnn 20 20 • Another remarkable feature is the behavior 
at small n, particularly at n = 0. Fig. 0] demonstrates 
that, e.g., the magnitude of Fq C 



000, 2m, 2m 



for m = 1 , 2 is 
only a few percent or even less of that of Gooo,2m,2m- 
Fig-Elshows that the ratio Fq 0u 21 21 /Goo n ,2i,2i becomes 
very small as n is lowered. A similar behavior has been 
found for all values (Xo,tp) on the glass transition lines 
we have investigated. The dips in F^ yz xx , /G xyz .xx' at 
n = demonstrate that the relaxation of the "self" part 
of the orientational correlators is practically not arrested 
by an orientational cage. 

Moving for oblate ellipsoids along the glass transition 
line towards the spherical limit Xq = 1, no qualitatively 
new behavior of the critical NEP is found, but it resem- 
bles always the characteristics of Figs. |3] and El However, 
this picture will change as we turn for oblate ellipsoids 
into the glass phase, as will be seen in subsection IIII CI 
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FIG. 7: Same as Fig. El but for prolate ellipsoids with a 
1.4524 and b = 0.72, i.e. (X ,ip) = (2.02, 0.394). 



I = V = 2, they are maximum exclusively at the zone 
center. A similar behavior is found for I ^ V (cf. Fig.^jJ, 
but here, e.g. for m = m! = 2, minima appear instead 
of maxima. None of the maxima of the static structure 
factors S l 2m,2m(q) and maxima/minima of 5 , 2m,4m(q) at 
the zone boundary persits in the limit t — > oo. Since 
these maxima belong to alternating orientational density 
fluctuations, this proves that such alternating local ar- 
rangements of the particles do not arrest. This can also 
be seen in real space. Fig. 01 exhibits the static orienta- 
tional density correlators G xyz ^ m ,2m and the F^ yz 2m 2m 
for a = 0.08 and b — 1.412. Indeed, the oscillations in the 
correlators Goon,20,20 and Go nn ,20,20 vanish completely in 
the long time limit, while the monotonous decay with n 
of the (m = ml > 0)-quantities is rather stable, even 
for infinite time. The almost vanishing of some critical 
NEP, however, does not require oscillations in the cor- 



The q-dependence of the critical NEP for prolate ellip- 
soids is sensitive on the shape of the ellipsoids (see Figs. 
[7| and ITUI) . The reader should note the higher values 
of the maxima in S\\i(q), which are necessary to get 
a glass transition, compared to the corresponding cor- 
relators for oblate ellipsoids. Let us have a closer look 
at prolate ellipsoids with a = 1.4524 and b = 0.72, i.e. 
(X ,cp) = (2.02,0.394). Fig.|7|shows that the structural 
arrest of these ellipsoids is completely different from that 
of oblate ones. The huge peak in ^2222(1 = (0,0, it)) 
at the zone boundary has a height of 157 and dominates 
the transition. Since this peak belongs to a wavelength 
of period two, it leads to strong frozen oscillations in the 
orientational density fluctuations on the lattice, as can be 
seen in direct space from Fig. |H1 Note that for the corre- 
lators Fq 0ti 21 21 almost no decay exists if n is increased. 
Again, like for oblate ellipsoids, the frozen F£ yz 20 2 o seem 
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FIG. 8: Same as Fig. 2] but for prolate ellipsoids with a = 
1.4524 and b = 0.72, i.e. {X ,ip) (2.02, 0.394). 



to play a special role, since they are much weaker than 
the NEP for m — m' > 0. Fig. shows the diagonal 
correlators for I = I' = 4. Note the very small scale 
for the static structure factors and NEP, in comparison 
with Fig. [7| Fig. [5] shows other interesting features of 
the MCT results for molecular crystals: besides the ap- 
pearance of simultaneous maxima of the normalized NEP 
at the zone center and its boundary (see /42,42(q) along 
the fourfold reciprocal space direction), the rule that the 
normalized NEP in reciprocal space are in phase with the 
corresponding static correlators |44| is violated. 

Finally, it must be said that the static structure fac- 
tors for ellipsoids with a = 1.4524, b = 0.72 have been 
calculated by OZ/PY theory. But MC results for 
other values of (a, b) in the vicinity of these parameters 
show that OZ /PY overestimates the maxima at the zone 
bondary in this region of the phase diagram. Therefore, 
the interpretation of Figs. EE] should be taken with some 
caution. Perhaps this overestimation is the indirect cause 
for the dip in ^ PY (Y ) for 2 < Y < 4 (see Fig.0). Why 
OZ/PY fails in this region of ellipsoids is currently un- 
known. 

As we turn to very elongated prolate ellipsoids, the 
transition scenario becomes simpler again. Fig. ED f° r 
a = 5.857 and b = 0.24 [yielding (Y , ip) = (24.4, 0.177)] 
serve as an illustration. The behavior of the (I = V = 2)- 
NEP with peaks at the zone center reminds one of the 
NEP for flat oblate ellipsoids (see Fig. 0). This means 
that for long prolate ellipsoids only nematic-like orienta- 
tional fluctuations may freeze. Such an extreme narrow- 
ness of the peaks at q = as seen in Fig. ^]is observed 
for prolate ellipsoids with Xq > 8 only, indicating the 
huge spatial extension of the frozen nematic-like fluctua- 
tions. 
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FIG. 9: Same as Fig. |S| but for prolate ellipsoids with a — 
1.4524 and b = 0.72, i.e. (X ,cp) (2.02,0.394), and for 
I = V =4, m = m' = 0, 1,2,3,4. 



C. Nonergodicity parameters in the glass phase 

In this subsection, we show by means of two examples 
how the NEP change in comparison to the critical NEP 
on moving slightly into the glass phase. The correspond- 
ing pairs (Xq, <p) are indicated in Fig. [21 too. 

For densely packed oblate ellipsoids with a = 0.78 and 
b = 1.1, i.e. (X ,ip) ^ (0.709,0.494) in the glass phase, 
the prototypical behavior of the critical NEP for oblate 
ellipsoids on the glass transition line shown in Figs. |3| 
and is clearly changed, as can be seen from Fig. ^] 
jS3j . Now, the Gaussian-like shape of the normalized 
NEP is much broader, indicating an enhanced arrest of 
orientational density fluctuations for q ^ 0, which is ex- 
pected due to the high packing fraction. This leads to 
a deviation of the frozen orientational correlators in di- 
rect space from the exclusive monotonous decay, which is 
present almost everywhere along the glass transition line 
for oblate ellipsoids. For example, the frozen i^oon 20,20 
for the (a, 6)-pair of Fig. ^2 (not shown here) have weak 
oscillations, reminiscent of the strong oscillations being 
present in the associated static Goon,20.20- 

Considering prolate ellipsoids with a = 1.7 and b = 
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0.66, i.e. (X 07 ip) ~ (2.58,0.389), slightly above the glass 
transition line, many different patterns of behavior occur 
in the NEP, as can be seen from Fig. This figure 

can directly be compared with Fig. [7J since the ellipsoids 
for both figures have almost the same packing fraction. 
Again, for one and the same NEP there partly exist si- 
multaneous maxima at the zone center and its bound- 
ary. Accordingly, in the limit of long times, we have 
frozen density-density correlators with either oscillatory 
or monotonous behavior, depending on AA'. 
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FIG. 11: Same as Fig. 01 but for oblate ellipsoids with a = 0.78 
and b = 1.1, i.e. (X ,<p) = (0.709, 0.494), above the glass line. 



D. Glass transition of the 1,1' odd correlators 

So far we have discussed the NEP for I, I' even. For I, V 
odd only the "self" part of the NEP is nonzero. It is use- 
ful to investigate the normalized, rotationally invariant 
"self" part of the NEP, i.e. 

As) _ Em=-I ^000, lm,lm 

■n — n ■ y m > 

/ 'm——t ^000, lm,lm 

Values for are given in Table [I] for those pairs (a, b) 
for which the glass transition has been found for I, I' odd. 
For comparison, the fj s ^ for I even are given in Table [IJ 
too. The relation 

// 8) </W, 1>V, (27) 
which is similar to 

f c (q) < f c (q') , Q>q' (28) 

for simple liquids, seems to be fulfilled for even and odd 
I separately. Note that the pairs (a, b) in Table |U are 
located in the glass phase for I, I' even. 



(a, fa) 


I = 1 


/ = 2 


I = 3 


I = 4 


(0.776, 1.1) 


0.203 


0.208 


0.124 


0.106 


(0.778,1.1) 


0.406 


0.268 


0.262 


0.143 


(0.78, 1.1) 


0.546 


0.333 


0.371 


0.185 


(1.7,0.66) 


3.00 x 10" 


2 0.197 


8.14 x 10 -3 


3.93 x 10~ 2 



TABLE I: Selected normalized NEP ft s > of the "self" part 
of the orientational density-density correlation function (see 
Eqs. (HJ). 
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IV. DISCUSSION AND CONCLUSIONS 

We have extended the mode coupling theory (MCT) 
for liquids to molecular crystals. The natural choice is 
the use of tensorial correlators, instead of correlators de- 
fined in a site-site representation j4f| . This leads for the 
dynamical correlators S\\> (q, t) to an integro-differential 
equation of second order in time. Truncating I at Z m ax, 
this set of equations is equivalent to the corresponding 
equation for a multi-component liquid of isotropic parti- 
cles (for binary systems see, e.g., |4(|). The memory ker- 
nel is approximated in the framework of MCT. The main 
differences to liquids are (i) the occurrence of Umklapp 
processes, if the sum of the wave vectors qi, q2 of the 
orientational density modes 5p\ x (qi) and 5p\ 2 (q 2 ) is out- 
side of the first Brillouin zone, (ii) besides the static two- 
point orientational correlators S\y(q) the need of the 
one-molecule orientational density pW(f2) as an input for 
the vertices of the memory kernel and (iii) the anisotropy 
of the static orientational current density correlators Jw 
which do not cancel completely from the memory kernel 
m\\'(c[,t). Nevertheless, the factor ksT/I of Jw drops 
out. Accordingly, the glassy dynamics and the ideal glass 
transition does not exhibit inertia effects, i.e. they are 
independent on I, the moment of inertia. Additionally, 
for rigid lattices, all Z,Z' = tensorial correlators vanish 
and can be skipped, due to the lack of TDOF. 

In order to discuss this set of MCT-equations for 
molecular crystals we have chosen hard ellipsoids of revo- 
lution with aspect ratio Xq = a/b fixed with their centers 
of mass at the sites of a sc-lattice with lattice constant 
equal to one. Increasing the size of the ellipsoids, which 
is equivalent to a decrease of the lattice constant, results 
in an increase of steric hindrance and finally in an ori- 
entational glass transition at the MCT-glass transition 
line (p c (Xo) shown in Fig. [2] for oblate and prolate ellip- 
soids. Since this orientational glass transition is mainly 
driven by the growth of S\\> (q) at the zone center or/and 
the zone edge, its origin lies in the growth of the orien- 
tational order close to but below the equilibrium phase 
transition line from OZ/PY theory. This is quite similar 
to what has been found for a liquid of hard ellipsoids [39j 
if the aspect ratio becomes larger than about 2 or smaller 
than about i. However, there is a difference between the 
molecular liquid (of ellipsoids) and the molecular crystal. 
Whereas the former already undergoes a glass transition 
for Xq > 2 or Xq < \ when 52m, 2m (0) is of order one, 
it must be S2m,2m(0) of order 10 for oblate (cf. Fig. EJ) 
or even order 100 for prolate ellipsoids (cf. Fig. 01. This 
proves that the translational degrees of freedom of the 
liquid still have a strong influence on the glass forma- 
tion, although they are not primarily responsible for the 
transition for Xq > 2 and Xq < |. This finding is consis- 
tent with results found from a MD simulation for difluo- 
rotetrachloroethane in its supercooled liquid and plastic 
crystal phase [43. For both phases, the glass transition 
temperatures T ( l iquid and T c plastic cr y stal were determined. 
That T c lk i uid = 139K > TP lastic CT ystai ^ 129K implies 



that the translational degrees of freedom of a liquid en- 
hance the glass formation which might be related to a 
facilitated cage formation for systems where the center 
of mass of the particles can move freely. 

Comparing S\\>(q) on the glass transition line for 
oblate (Figs. |3| and EJ) with prolate ellipsoids (Figs. 
and I9I1U[I already shows that the tendency to an ori- 
entational glass formation for oblate ellipsoids is larger 
than for prolate ones. This can also be seen from Fig. [21 
since the distance ippy(Xq) — ip c (Xo) is large for very flat 
oblate ellipsoids, only. This difference may be explained 
as follows. If we fix the length a of prolate ellipsoids 
and decrease their thickness b to zero, then the excluded 
volume interaction becomes zero. Particularly, the static 
correlators become trivial, which leads to vanishing ver- 
tices and consequently to a disappearance of the MCT 
glass transition. If on the other hand we fix the diame- 
ter b of the oblate ellipsoids and decrease their thickness 
a to zero the excluded volume interaction still exists for 
a = 0. This seems to be an important difference between 
oblate and prolate particles. 

From the solution of the MCT equations for t — » 00 
we obtained results for the critical NEP F xx , (q) and the 
corresponding normalized ones, /£y(q), as well as NEP 
deeper in the glass. Due to the lattice translational in- 
variance, q can be restricted to the first Brillouin zone. 
Within this zone the critical NEP do not have much 
structure. Almost all of them either exhibit a maximum 
(or minimum for I ^ I') at the zone center and/or at 
its edge, depending on AA' and the direction of q. How- 
ever, going deeper into the glass and varying the ellipsoid 
shape, i.e. Xq and/or ip, leads to significant changes in 
the q-dependence especially of the normalized NEP, as 
demonstrated in Figs. Illlandll2l 

The MD-simulations for cyanoadamantane 0| and 
chloradamantane |18| reveal quite similar glassy dynam- 
ics as found for supercooled liquids Q, H, H • Particularly, 
the authors of Refs. 0, 0] stress that their molecu- 
lar crystals can be supercooled and that the relaxational 
dynamics is consistent with MCT predictions, at least 
where this has been checked 0, . Since both model 
systems exhibit tremendous slowing down in the super- 
cooled regime an orientational glass transition or at least 
a crossover from ergodic to quasi-nonergodic dynamics 
must also occur in the supercooled phase. This is dif- 
ferent from what we have found for the hard ellipsoids. 
In our case the glass transition line is located within the 
dynamically disordered equilibrium phase of PY theory, 
which itself is a consequence of the static input taken 
from PY theory. It would be very interesting to insert 
the static correlators [5^ from the MD simulations in the 
supercooled phase into our MCT equations and to check 
whether one obtains a glass transition. The investigation 
of hard ellipsoids on the sc lattice has demonstrated that 
the magnitude of the extrema in S\y (q) at the zone cen- 
ter or/and edge must be rather large (at least for prolate 
particles). It is not obvious that the simulational results 
in the supercooled phase fulfill this criterion. Of course, 
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it could be that the average height of S\y (q) in (/-space is 
much larger than for the ellipsoids such that large max- 
ima/minima of S\\i(q) at the zone center and/or edge 
are not really necessary. We hope that these questions 
can be answered in future. 

To conclude, we have shown that MCT can be derived 
for molecular crystals, as well. Whether or not the MCT 
approximation (which mainly consists of the factorization 
of a time dependent four point correlator) is also a rea- 
sonable approximation for molecular crystals as it is for 
glass-forming liquids has to be investigated by compari- 
son of the results from present MCT for molecular crys- 
tals with simulational and experimental results. As al- 
ready mentioned above our conventional MCT approach 
will become worse if the thickness of prolate particles be- 
comes small. In that case it is the entanglement which 
is responsible for glassy dynamics [22T 123^ . This requires 
a different theoretical description, as recently discussed 

MM. 



Since ui n ■ L^„ = <-^' n ■ L^ (where the primed quantities 
refer to the body fixed frame) we get for the canonical 
average in Eq. IjAljl in close analogy to molecular liquids 

M 



oca' 

aa' 

(A2) 



where, a = x,y, z are the cartesian components in the 
body fixed frame. This leads to 



APPENDIX A: CALCULATION OF J AV (q) = J\\> 

Substituting the A-Fourier transform j\(q) of j n (Q,) 
(Eq. 10) into Eq. yields 



kT 1 

J xx ,(a) = Aii — i l - l ( (La n Y x {Sl n ))* ■ {-Ln n Y y {n n ))) , 

(A3) 



JAv(q)-^^^E eiq ' x ""' x 

„ » which is q- and n- independent, and can therefore be eval- 

x ( (u n ■ Lfi n Y A (fi„))* (uv • L Qn ,Yy(fl n ,)) ) . (Al) uated for ar bit ra ry n. 



Using L^Yimip*) — ^ m i = _iLf mm ,Yi m i{VL) and the product rule for the spherical harmonics and substituting the 
explicit expression for Lf mm , we get with c(l, I', I") as in ljC4() 



J AA /(q) = 4tt — i 1 '- 1 (-l) m [ mm ' -mm'm") 



-yjl(l + 1) - m(m + 1) + 1) - m'(m' + 1) C(llT, -m -l,m'+ 1, m") 



-y/l(l + 1) - m(m - 1) y/V{V + 1) - m'(m' - 1) C(llT, -m + 1, rri - 1, m") J c(ll'l") ( Y\» ) . (A4) 



This expression strongly simplifies since 



1 



[• • • ] Clll'l") = -- V / W+l) y/l'(l' + l) 



This leads to the final result 
kT 



(2Z + 1)(2/' + 1) 



4tt (21" + 1) 



{CQl'l", 1 - 10) + C(ll'l", -110)} C(llT, -mm'm 1 ' 



(A5) 



JAA'(q) =^ — l l - l (-l) m+1 - y/l(r+rjy/l'(l' + l) 

"(2/ + 1)(2Z' + 1) 



E 



4tt (21" + 1) 



{CQl'l", 1 - 10) + C(ll'l", -110)} C(ll'l", -mm'm") ( Y\„ ) . (A6) 
I 



Note that (Y\) is given by 



(y x )= I <m P W(n)Y x (n) = (-i) l p[ 1) , (a?) 



i.e. J\A'(l) — ^aa' only involves the A-transform of 
pM(Q). 
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APPENDIX B: MODE COUPLING 
APPROXIMATION 

In this appendix we shortly describe the mode cou- 
pling approximation leading to the results presented by 
Eqs. Cljl-lQSJ). 

The derivation of the Mori-Zwanzig equation is stan- 
dard by using the projectors onto the slow variables 
5p x (q) and j A (q), 

A-7T 

P e = 7fY,'\ ) ( S_1 (q))^ < ^v(q) I . (Bl) 



4-7T 

^ = ^E1^)>( J ~V'(iA'(q) 



XX' 



N 



(B2) 



XX' 



The prime on sums denotes summation such that I, I' > 0. 
The projector Q in Eq. Q14[) is the given by 



Q = l-P p -Pj 



(B3) 



In order to approximate M\y (q, t) Eq. I|14|) we introduce 
the projector onto pairs of orientational density modes: 

V = E E' |^A 1 (qi)^A 1 (q' 1 )>x 

qiqi,q2q 2 AiAiA 2 A 2 
e i.bz 

x .9A 1 A' 1 ;A 2 Ai(qiqi,q2q , 2 ) <*PA,(<te)*PAj(qa) I . ( B4 ) 

where (<7AiAi,A 2 A 2 (qiqii q2q^)) is the inverse of the static 
four-point correlation matrix (( 6p x (qi) 6p x , (q^) x 
x <!)p A2 (q 2 ) #Pa' (q'2) ))• We use the approximation 



.9A 1 A' 1 ;A 2 A^(qiqi,q2q 2 ) 



1 /4tt 
4 ViV 



(B5) 



<W <5 q ;,q 2 (S- 1 (qi))A 1 A 2 (S- 1 (q 2 ))A i A i + (1 - 2) 



which is consistent with the mode coupling approxima- 
tion of Mx\'(q,t) for t = (see Eq. (JbTJ). 

The mode coupling approximation consists of two main 
steps. First, the fluctuating force is approximated 



Q I £jA(q) ) ~ VQ I Cj\(q)) 



(B6) 



Substituting (|B6I) into Eq. I|14l) leads to a time-dependent 
four-point correlator, which in a second approximation is 



factorized. For qi,q' 1 ,q2,q2 £ l.BZ we have 

< 5p* Xl (qr) Sp* x , (q[)\Qe- lQCQt Q I <5p A2 (q 2 ) 5p A , (q' 2 ) ) 



N 2 

With these approximations we obtain 



Sq uqa Vi.qi SAiA 2 (qi>*) ^a;^ (q2,*) + (1 2) 

(B7) 



MAA'(q, t) 



1 /4tt 



E E' 

eiSz x } x f x f x f 

X X X 2 X 3 X 4 



x ((^jA(q))* |Q|<5p Al (qi)£p A2 (q 2 )) x 

x (S- 1 (q 1 )) AlA3 (S- 1 (q 2 )) A2A4 x 

x S A3 ^(qi,t) S^fat) x 

x (S- 1 (q 1 )) A , Ai (S- 1 (q 2 )) A , A ,x 

x <^ A ,(q 2 )5p A ,(q 1 )|Q|£iv(q)>. (B8) 



APPENDIX C: CALCULATION OF 

((A7A(q))* IQ!5pA 1 (qi)<5pA 2 (q2)} 

This correlator is calculated quite similar to simple [j| 
and molecular liquids [28|, |23, H3 by using Eq. (|B3(I and 
Pj I Sp\ x (qi) <5/9 Aa (q 2 ) ) = 0, due to time reversal symme- 
try. Then we get 

( OOA(q))* I Q I Wqi) Wq 2 ) ) = 

= <(£j A (q))* ( 5p Al (q 1 )^ A2 (q 2 )) 

- < (Cjx(q)T I P P I Sp Xl (qi) ^A 2 (q 2 ) ) • (CI) 

The first term on the r.h.s. of Eq. (|C1|) is easily rewrit- 
ten by taking into account the hermiticity of C and the 
continuity equation Eq. (0} and Eq. ©. This leads to 

((47A(q))*£pAi(qi)5pA a (q2)) = 

= ( Jjt(q) Jx 1 (qi) 5pa 2 (qa) > + (1 <- 2) . (C2) 
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Substituing the A- Fourier transform of 6p n (Q) and jn(f2) into the r.h.s of Eq. (|C2|) we arrive at 



N kT 



m-\-m 



Q A" 

-- \/Z(Z + 1) - TO ( m + 1) yjh{h + 1) - mi (mi + 1) C{ll-yl" ,-m - I, mi + 1,-m") 



mmi C{ll\l , —rami — m ) 



+ 1) - m(m - 1) + 1) - mi(mi - 1) Cfjhl", -m + l,m 1 - 1, -m") 



with 



c(ZZ'Z") 



(2Z+1)(2Z' + 1) 
4tt (2Z" + 1) 



C(llT, 000). 



c(lhl")Sy, X2 (<i 2 ) (C3) 



(C4) 



Here we have used the product rule for the spherical harmonics and the factorization of canonic integrals as in Eq. 
(|A2|) . The second term on the r.h.s. of Eq. (|C1|) is rewritten by using P p from Eq. IjBlfl and again the hermiticity of 
C, as well as the continuity equation: 



((£jA(q))* \P P \6pxd<ix)Sp)»{<to)) = E' (S- 1 (q))vA» <<5p A »(q)^A 1 (qi)^A 2 (q 2 )) 



(C5) 



A' A'' 



Substituting Eqs. (|C2|) . (|C3|) and (|C5|) into Eq. I|C1|I . the l.h.s of Eq. I|C1|1 is expressed by the static two-point and 
three-point correlators SxA'(q) an< ^ (^Px"( c l) <^°Ai(qi) ^Pa 2 (Q2) ), respectively, and by Jaa'- J = (^AA') is calculated 
in App. El ( 5p* x „ (q) Sp Xl (qi) ^a 2 (qa) ) in App. [Q1 
Now we rewrite J2\"' ■ ■ • in Eq. (|C3|I as follows: 



E'--- = E'E'* Zl 



(-i) 



m+m ' 



A" A', 



mm 1 C(ll 1 l ,—mm 1 —m ) 
+ 1) - m(m + 1) ^i(/i + 1) - mi (mi + 1) C(/Z^", -to - 1, mi + 1,-m") 
i V^G + 1) - ™(™ - 1) + 1) - mi (mi - 1) C(ZZiZ", -to + 1, m' a - 1,-m") 

E (S- 1 (qi))A' 1 A'" ^A'"A 1 (qi)5 A » A2 (q 2 ) , 



and substitute succesively the terms on the r.h.s. of 



S- 1 (qi)=d- 1 -d- 1 +D- 1 -— c(qi), 



(C6) 



(C7) 



which is a rearrangement of the OZ equation 1)20(1 . into Eq. i|C6(l . 

In the first step we replace d^ y „ Sy/\ 1 (qi) S\"x 2 (q 2 ) = c?a"a^ ( 5 )a; a x (qi) ( ) a ^a 2 (qO- This expression 
arises if the matrix elements of d 1 on the r.h.s. of 1C7JI are used with l|C6(l . Using the explicit result for the matrix 
d (see and the relations 

[mmi c(ZZ' x /") C(ll[l", -mm' x m n ) c(l"l' 2 l') C(l"l' 2 l', m"m' 2 m') 

A" 

+ toto' 2 c(ll' 2 l") C(ll' 2 l", -mm' 2 m") c(l"l[l') C(l"l[l', m"m' 1 m l ) ] 
= ^ mm"c(l[l' 2 l") C(l[l' 2 l", m! x m! 2 m") c(ll"l') C(ll"l', -mm"m') , (C8) 



E [y/l(l + l)-m(mTl) sJWi + 1) - mi(mi T 1) c(^Z") C(Z/iZ", -to ± 1, mi T 1, m") c{l"l' 2 l') C(l"l' 2 l', m"m' 2 m') 

A" 

+ + 1) - to(to t 1) \Jl' 2 ( l 2 + 1) - ™ 2 (™ 2 T 1) c(ZZ 2 C^", -m ± 1, "4 T 1, m") c(/"/i/') C(l"l' x l\ m"m[m') ] 

= ^ + 1) - m(m t 1) \/Z"(Z" + 1) - m"(m" T 1) c(ZiZ 2 Z") C(l[l' 2 l", m' x m' 2 m") c(ll"l') C(ll"l', -m ± 1, to" t 1, m') . 

(C9) 
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we find that this part of (|C6|I taken together with the same part in the partner expression of i|C6|) due to l|C2|l cancels 
with the part ( (£j A (q))* \ P p \ Sp Xl (qi) Sp X2 (q 2 ) ) of JOU, if Eqs. CSS and are used in Eq. 

We turn to the term D _1 - d 1 of (|C7) . which leads to (D _1 - d~ 1 )A< L A""S , A'"A 1 (qi) if substituted in l(C6)l . Since 



[D^-d- 1 ] (fl.fi') 



4tt 



4tt p( x )(fi) 4tt pW(fi') (4tt) 2 



p«(fi) 



dfi 



(CIO) 



D 1 d 1 consists just of a nontrivial first row and column, while S(qi) has a vanishing first row and column. So the 
product (D _1 — d _1 ) S(qi) has nonvanishing elements in its first row, only. Therefore, (D _1 — d -1 )^ A""Sv"Ai (q x ) = 0, 
if not l[ — m'i — 0. But if we evaluate the coefficients of l|C6ll with l\ — m[ = 0, we find that the part D 1 — d _1 of 
(|C7|) contributes nothing. 

What remains is the last term on the r.h.s. of l|C7(l . If substituted into Eq. (|C6|) and the corresponding partner 
expression due to Eq. I|C2J) . respectively, this term delivers the final result 

<(£jA(q))*|Q| Wqi)Wq 2 )> = 

A J" E<W<^+Q E ' [^20C(/^/,m! i m'2^)^Ai(qi)5 A;Al (q 1 )5 A , A2 (q 2 ) + (1 <- 2) ] , (Cll) 



Q 



with b(l,l',l") from Eq. (|18[) . Here we have used the relation l|A5(l for the Clebsch-Gordan-coefficients. If Eq. I|C11|) 
and its conjugate is substituted into Eq. (|B8|) one obtains the mode coupling approximation of the slow part of 
M AA /(q, t), which then leads to the final result for m(q, t), Eqs. (fH)l)-((T§)l. 



APPENDIX D: APPROXIMATION OF 

( 5p\ (qi)<5pA 2 (q2)5pA 3 (q3) } 



The approximation of the static three-point correlator 
is rather involved. Therefore, the most crucial steps are 
presented only. Readers which are interested in more 
details are referred to Ref. 01 • 

The corresponding static three-point-correlator for 
simple liquids was approximated by the convolution ap- 
proximation 5] . It has been proven that the approxima- 
tion of the corresponding correlator for molecular liquids 
|28j i s ag ain the convolution approximation as defined in 
Ref. |50| . However, performing the convolution approx- 
imation for molecular crystals does not lead to a simple 
result. Therefore, we have chosen a different approx- 
imation. ( Sp* Xi (qi) 5p\ 2 (q 2 ) 8p\ 3 (q 3 ) ) is the A-Fourier 
transform of ( Sp ni (fii) <5p„ 2 (fi 2 ) 5p n3 (fi 3 ) ) given by: 



<^ Al (qi)^A 2 (q2)<5pA 3 (q3)) = 



E 



'i(-qi-x„ 1 +q 2 -x„ 2 +q3-x rl3 ) -h+h-h x 



(Mi dn 2 dn 3 y A * (no y A2 (fi 2 ) y X3 (fi 3 ) x 



x (5p ni (fii)Sp n2 (n 2 )Sp na (fi3)) 



(Dl) 



For ( Sp ni (fii) 8p n2 (fi 2 ) 5p n3 (fi 3 ) ), one can prove that a 



reasonable approximation is 
(WOi) WOa) -W^)) « E / dn P (1) (n) 

n J 

, £7) G„„ 2 (fi, fi 2 ) G nn3 (Vl,Vl 3 ) 



p(D(n) p«(fi) pW(o) 



(D2) 



where G„ ni (fi, fii) = G ni „(fii,fi) has been used. Per- 
forming the Fourier sums of Eq. i|Dl|) on approxima- 
tion lf02"|l yields 



N 



q2+q3,qi+Q 



Q 



dfip (1) (fi) x 



s(gi,fii,fi) s(g 2 ,fi,fi 2 ) s(g 3 ,fi,fi 3 ) m ^ 
p(!)(Q) p(D(n) p(!)(fi) ' 1 j 



where 



5(q, fi, fi') = 4tt ^ G nn , (fi, fi') e ,qx -' . (D4) 

x nn / 

Substituting 

5(q,^fi')=E'H) i '^AA'(q)^(0)F A *,(n / ) (D5) 



AA' 



and 



- 4tt J dn"d-\n,Q")S(< l ,n"M) = 
= 47r E E '(-*)''"' aa' (<0 *M >a* . 

A A' 

(D6) 
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with 



4vr pW(fi) 



(D7a) 
(D7b) 



into Eq. l|D3Jl and taking the A-transforms as defined in 
Eq. I|D1I) afterwards we get 

( s p*x 1 (qi) s px 2 (q2) Spx a (q 3 ) > « 



4tt 



Q 



Although, the product of the last three factors of 
Eq. (|DB|> does not look symmetric, one can show that 
all three factors indeed are equivalent. 



x c(l 2 l 3 l 1 )C(l 2 l 3 l 1 ,m 2 m 3 m 1 )x 



S Al Ai(qi) (d- 1 5) A , A2 (q 2 ) (^ 1 5) A , A3 (q 3 ). (D8) 
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